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The density profile and surface tension for the interface of phase-separated colloid-polymer mix- 
tures have been studied in the framework of the square gradient approximation for both ideal and 
interacting polymers in good solvent. The calculations show that in the presence of polymer-polymer 
excluded volume interactions the interfaces have lower widths and surface tensions compared to the 
case of ideal polymers. These results are a direct consequence of the shorter range and smaller depth 
of the depletion potential between colloidal particles induced by interacting polymers. 

o . 

O . I. INTRODUCTION 

Oh- Colloids are "soft" materials, readily deformable, with weak interfaces. This can be easily inferred from the "giant 
^ i atom" picture of colloidal suspensions where, even though each colloid is made up of thousands of individual molecules, 
it is treated as a single particle interacting with an effective potential^. Since the effective interactions are roughly 
of the same shape as those of atomic fluids, an approximate corresponding states principle suggests that the reduced 
or dimensionless surface tensions should be similar. Near the fluid- fluid transition, the attractive interactions for 
both classes of fluids are typically of order fceT, but the colloidal particles have radii R c which can be 2 or 3 orders 
of magnitude larger than molecules. Thus the surface tension, which scales as 7 ~ fc^T '/ 'R%, is expected to be 4 or 
more orders of magnitude lower than the values found for simpler atomic and molecular fluids. Similar approximate 
corresponding states arguments also explain why colloidal crystals are so easily deformable: their elastic constants, 
which scale as ksT/R^., are at least 6 orders of magnitude lower than those of simple atomic or molecular crystals. 
Colloids are indeed a form of "soft matter" . 

Surface tension plays an important role in the formation of interfaces, as well as in phase transition kinetics, 
1-^ ' nucleation and spinodal decomposition^. Its indirect effects are therefore easily observable, but its low values make 
direct experimental measurements very difficult. Nevertheless, some recent experiments have made significant progress 
in measuring the fluid-fluid interface of colloid-polymer mixtures and its surface tension^. In these systems, adding 
O , non-adsorbing polymers induces attractive depletion pair potentials between the colloids^, which lead to the observed 
phase-transition between a colloid-rich ("liquid") and a colloid-poor ("gas") phase, separated by an interface. Because 
t—I [ the experimental parameters can be easily tuned and controlled, colloid-polymer mixtures form an important model 
■ system for the study of phase transitions in soft matteA 

By applying theories similar to those used successfully for atomic and molecular fluidsiS, Vrijii, and Brader and 
Evans^, calculated the properties of this fluid-fluid interface for the case of ideal polymers, finding qualitative agree- 
ment with experiments. We have recently derived a depletion pair potential valid for interacting polymers^, which 
captures the dominant effects of polymer-polymer interactions on the phase diagrams^. This success suggests that 
• the same potential can be used to calculate the properties of the fluid- fluid interface. 

The main purpose of this paper is to investigate the effects of polymer-polymer interactions on the fluid-fluid 
interface of colloid-polymer mixtures. For that reason, we apply the same combination of thermodynamic perturbation 
theorjii^ and square-gradient theory^ that was used by Brader and Evans^, but with the new potential^ instead of 
the Asakura-Oosawa (AO)£i£ pair potential, valid only for ideal polymers. The differences between our new results, 
' and those of ref£, are then mainly due to the effect of polymer-polymer interactions. 

The use of colloid-colloid depletion pair potentials describes one level of coarse-graining. It is also possible to 
derive a more fundamental two-component picture based on polymer-polymer, polymer-colloid, and colloid-colloid 
pair potentials. A number of more recent investigations have used sophisticated two-component density functional 
(DFT) theories for the AO modeli^ to uncover a host of interesting interfacial phenomena, including oscillatory density 
profiles at the fluid-fluid interface and a series of layering transitions at the fluid-hard-wall interfacei&. Computer 
simulations^ have confirmed some of these results. At present, all these theories are only applicable to the AO model, 
and it is unfortunately not yet clear how to extend them to interacting polymers (see however ref^) For that reason 
we restrict ourselves to the simplest square gradient approximation for the interfacial profiles. 

Our paper is organised as follows: after briefly reviewing the nature of the depletion potentials and the equilibrium 
phase-diagrams in section II, we describe the implementation of the square gradient approximation in section III, and 
present our results for the interfacial tension and width in section IV. Section V summarises our conclusions. 
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II. ONE-COMPONENT EFFECTIVE DEPLETION INTERACTIONS 



This section briefly describes the colloid-colloid effective depletion potentials for both ideal and interacting polymers. 
They are characterised by the polymer radius of gyration R g , the colloid radius R c , and the polymer number density 
p p , or equivalently by the size-ratio q = R g /R c and the reduced polymer density r/ p = p P /p P , where p* — ^R g is the 
so-called overlap density. In the ideal case, the depletion interaction between two isolated colloidal spheres at distance 
r is accurately approximated by a potential of the Asakura-Oosawa form: 
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for 2R C < r < 2{R C + R e J^); V id (r) = for r > 2(R C + R e J£). Here, a cp = (R c + R e J£) and p r p is the polymer density 
in a reservoir in osmotic equilibrium with the full colloid-polymer mixture 19 . The range of this potential depends only 
on the polymer length and the depth is proportional to the polymer density. The effective Asakura-Oosawa parameter 
R e /o is set by the requirement that the insertion free energy of one colloid is equal to that of ideal polymersSS; it is 
given by2i 
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For interacting polymers, we will use a recently proposed pair potential^, which accurately reproduces the depletion 
potentials obtained from direct computer simulations: 



V s (r) = -nR c7w (p r p )D s ( P ;) 1 - 
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for 2R C < r < 2R C + D s and V s (r) = for r > 2R C + D s . Here 7ui(pp is the surface tension of the polymer solution 
near a single wall^i and D s {p r p ) is the range of the potential, given by 
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where n is the osmotic pressure of the solution of interacting polymers, which is well understood^. 

The range of Va is independent of density, whereas the range of V s shrinks with increasing density. Furthermore, 
for a given p r p and R g , the well-depth of Vid is greater than of V s , which implies that ideal polymers induce stronger 
depletion potentials than interacting polymers (see e.g. Fig. 2 of ref(i£ for some explicit examples). 

The differences in pair potentials help explain why, for a given q. phase-separation occurs at a larger value of rj p 
for interacting polymers than for ideal polymersSi^Si, a difference that grows with increasing q. Of course the pair- 
potential approximation becomes increasingly unreliable at high q values, as many-body interactions become more 
important. However, we have recently showni^ that calculations based on pair potentials alone remain remarkably 
accurate up to q ~ 1. In this paper we use the phase-diagrams calculated in ref^, based on second order perturbation 
theory, as the basis for our calculations of the properties of the fluid-fluid interface. We make the implicit assumption 
that the effective Hamiltonian used for phase-behaviour is also appropriate for describing the interface. This follows 
from the fact that we are working at a contact chemical potential, so that the same effective potential holds across 
the density inhomogeneity occurring at the "free" interface (see ref^). 



III. INTERFACIAL PROPERTIES FROM THE SQUARE GRADIENT APPROXIMATION 

Once phase separation occurs, there are two phases with well defined colloidal bulk densities (pf and p^ for the 
dilute and concentrated colloidal phases, respectively). Both phases are separated by a planar interface where the 
local density depends on the distance to the interface, p c (z). A well known way to treat the free-energy cost of making 
an interface is given by the square gradient approximation, where the free-energy is expanded to lowest non-trivial 
order in a gradient expansion around the homogeneous fluid. The surface tension and the density profile are then 
obtained from the integral of the free-energy across the interfacei2i2S*2LS£ 

dz, (5) 
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where '3/ (p c (z)) = f (p c (z)) — p c Pc(z) + F '■ Here p, c and P are the chemical potential and osmotic pressure of the 
colloids at coexistence and / (p c (z)) is the Helmholtz free energy density of a hypothetical colloid fluid of density p c (z). 
The coefficient of the square gradient term, n, describes the free-energy penalty for creating an interface. Minimising 
this functionalist leads to the following expressions for the density profile 

az J k 
and the surface tension 

7 = 2/ [n^] 1/2 dp c . (7) 

Requiring the functional in Eq. JSJ to satisfy linear response relates the coefficient n to properties of the direct 
correlation function c(r) of the homogeneous fluicUS 



r 4 c(r, p c )dr. (8) 



Note that all these variables depend implicitly on the polymer chemical potential (or equivalently the polymer 
reservoir density, p r p ) of the corresponding coexistence point. Due to the factor r 4 in the integrand of expression JSJl, 
the value of k is mainly determined by the behaviour of c(r) at large r, where it is well known that c(r) rj —j3V{r). 
We therefore follow Refi^, and set c(r, p c ) to be zero for r < 2R C and equal to —(3V(r) for r > 2R C , where V is given 
by Eqs. 0) and © for ideal and interacting polymers, respectively. Hence Eq. JHJ reduces to 
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This approximation has the further advantage that it circumvents the conceptual difficulty of defining c(r; p c ) in 
the coexistence region. Even though the approximation for c(r) itself may not always be so reliable, we found that 
the values of k still compare well with more sophisticated calculations of c(r), because this simple model interpolates 
between the values at the two coexistence points. Similar conclusions were reached in a paper studying Lennard Jones 
systems^, where the Percus-Yevick approximation for the low and high density fluid phases was combined with a 
lever rule to obtain c(r) in the coexistence region. 



IV. RESULTS AND DISCUSSION 



A. Phase behaviour and interfacial properties 

In our earlier work 1 ^, the free energy densities for the effective one-component system, f(p c ), were calculated for 
various q-ratios by second order perturbation theory using the Barker-Henderson formulation 31 . The phase diagrams 
were determined by the common tangent construction. The resulting coexistence curves for ideal and interacting 
polymers are plotted in Fig. for size-ratio q = 0.67, as a function of the colloid packing fraction r\ c = Airp c R^/3 
and polymer reservoir packing fraction rf v = AirppRg/3. The fluid- fluid coexistence lines are at higher values of 
polymer packing fraction for interacting polymers, which implies that polymer excluded volume effects reduce the 
global attraction between colloidal particles. But not only is the position of the binodal different, so is its shape. In 
particular, the binodal is flatter and the separation between the critical and triple points is smaller for interacting 
polymers. This effect is not merely an artifact of the pair-potential approximation, since the qualitative difference 
in shape is also observed when comparing two-component simulations of ideal and interacting polymer models^. 
Although we only show results for one size-ratio in Fig. ^ the differences become more pronounced for increasing 
qii^iSi, and are finally quite dramatic in the so called "protein limit" where q >> 1—. 

It should be noted that, at least within our perturbation theory treatment, the gas-liquid binodal obtained for 
interacting polymers at q = 0.34 is metastable with respect to the fluid-solid coexistence. This is not, however, an 
obstacle to the calculation of surface tensions. Furthermore, in many experimental systems the fluid-solid nucleation 
rates are very low, allowing the observation of metastable fluid-fluid phase-separation. 

The free energy densities from perturbation theory were used in Eqs. © and Q to calculate first the density profiles 
and then the surface tensions for the coexistence points along the fluid-fluid binodal. Two typical density profiles 
corresponding to ideal and interacting polymers are shown in Fig.|21for q = 1.05 and Ar^ = (rf — vjl' cr%t ) / 'tf ,c "* = 0.2. 



All profiles obtained using the square gradient theory share approximately the same shape, i.e. a smooth monotonic 
curve which goes from the dense to the dilute colloidal phase. Comparison between both curves shows that the 
difference between the colloidal packing fractions in the two phases is larger for interacting polymers (a consequence 
of the flatter binodals), whereas the interfacial thickness is smaller. 

The interfacial profiles are essentially characterised by their width. The 10-90 width of the interface (W), defined as 
the distance along the interface over which the colloidal density varies from (77^ + 0.1(77^ — 77^)) to (77^ + 0.9(77^ —77^)), 
is plotted in Fig. I^w.r.t. the deviation of the polymer reservoir packing fraction from the critical point. As expected, 
the interface is very diffuse near the critical point, and becomes sharper upon approaching the triple point. The 
interfacial widths increase with increasing q, reflecting the longer ranged attractions. For a given value of (7, the 
widths are consistently lower for interacting polymers than for ideal polymers. 

The resulting dimensionless surface tensions (jR^/kBT) are shown in Fig.0]versus the difference in colloidal packing 
fractions between both phases for three different size ratios, q = 0.34, 0.67 and 1.05. For both ideal and interacting 
polymers, the surface tension vanishes at the critical point, and increases with At7 c = 77^ — 77^ as expected. These 
reduced values are of the same order as those found for simple liquids^. Their absolute values depend only on R c ; 
typical colloid sizes result in values of 7 near the triple point on the order of jj,N/m which is much smaller than the 
values of 10 — 100mN / m found for simple fluids. For a given A?7 C , the surface tension increases with q. Moreover, for 
a given q the ideal polymer surface tension is always significantly larger than that of interacting polymers. 

In view of the mean field nature of the present theory, the critical exponents are obviously classical, i.e. 7 cx Arj^ 
and W tx A77" 1 . However, the prefactors are expected to be different for ideal and interacting polymer depletants. 
These prefactors are determined in the following sub-section. For ideal polymers similar scaling laws with the reservoir 
density p r p can also be derived, by exploiting the analogy with inverse temperature. For example^, W cx (Arjp)~2. 
However, for interacting polymers such a simple scaling with p r p does not follow, because the effective well-depths (i.e. 
the inverse "temperatures") don't scale in a simple way with this variable. 



B. Connection to the form of the depletion potential 

The results of the previous subsection show that both the the 10-90 width, and the interfacial tension, are lower in 
the case of interacting polymers. In this section, we will attempt to rationalise these differences on the basis of the 
effective pair potentials. 

Close to the critical point, ^(rj c ) can be approximated as the product of two quadratic potential wells centred 
around the coexistence points^ 

^ (r)c-V?) 2 (r]c-tf) 2 , (10) 

where \&* = (4ttRI/3)(3^> is dimensionless and C = (1 / 12) d 4 ^> * /dr] 4 , calculated at the critical point. If we assume 
that k is independent of p c , as we did in the previous section, then inserting Eq. 110(1 into Eq. J7J) yields the following 
expression for the reduced surface tension: 
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where k* = /3k/ R^.. The dependence of k* on the polymer reservoir packing fraction is indicated explicitly to remind 
us that it is not a constant parameter but rather increases as we move away from the critical point since it depends 
on 77p through V(r). The same arguments lead to£ 
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The values of W and 7* are determined by the parameters k* and C. Eq. J^J directly relates k to the depletion 
potential, whereas C is given by the shape of the free energy inside the van der Waals loop. Understanding how the 
depletion potentials govern the interfacial behaviour now reduces to explaining how these two parameters depend on 
the potentials. 

Fig. compares the depletion pair potential for ideal and interacting polymers at their corresponding critical points 
when q = 0.67. Even though the depletion potential for interacting polymers has a larger depth at contact, its range 
is significantly shortened by the polymer excluded volume interactions^. The integrand in Eq. which determines 
k, multiplies the potential by r 4 , giving extra weight to the longer ranged parts of the potentials. Thus the ratio of 
K* nt (calculated with Eq. ©) to K* d (calculated with the potential of Eq. Q), is always less than 1, and decreases 



with increasing q, as shown in the inset of Fig. [5] It has been recently shown that the reduced second virial coefficient, 
which is proportional to the integral of r 2 (exp[— f3V(r)] — 1), is very similar at the critical point for a wide class of 
attractive potentials^. This observation remains true for the two depletion potentials, giving further support to our 

argument that it is the large factor r 4 in the integrand for k which is responsible for the differences between n* d and 

* 

K inf 

The value of the other parameter, C, is not as easy to link directly to the pair potential. Nevertheless, it can be 
determined numerically from our previous calculations of the surface tension. It turns out not to depend strongly on 
q. We find C = 85.8 ±0.5 for ideal and C = 76.0 ±0.3 for interacting polymers, the small difference perhaps reflecting 
the fact that the free energy loop is slightly flatter for interacting polymers. 

Compared to the large changes in the surface tension, the values of C are quite similar, so that differences in the 
7*s, which scale as 7 ~ v k*C, arise mainly from k. Since k increases with the range of the potential, i.e. with q, 
this explains why, for a given type of depletant, the surface tension grows with q. It further shows how the main 
differences in 7* between ideal and interacting polymers are linked to the reduction of the depletion potential range 
induced by polymer-polymer interactions. Since this effect becomes more important for increasing q, we expect the 
differences between the surface tensions induced by ideal and interacting polymers to grow with q as well, and to 
become more pronounced in the protein limit^, as was recently pointed out by Sear^Si. 

The interfacial width scales as W ~ v/ k*/C. Again, differences in W are dominated by k. Thus, the net effect of 
adding polymer-polymer interactions is to decrease the interfacial width. 

These results show that the decrease of the range in the depletion potential, caused by the polymer interactions, 
plays the dominant role in determining the differences in the interfacial properties between the two types of depletants. 
The change in the depth of the potential is only a secondary effect. 

In all the arguments above it should be kept in mind that our double symmetric parabola approximation of Eq. (jl 1|) 
is only valid close to the critical point. For coexistence points far from it, the bulk correlation length (defined as 
£b = (2k/ (dfic(pb) I dpb)) 1 ^ 2 , where pb is the colloidal density in the bulk phase) does differ between the two phases, 
implying that the decay of the profile tails is different between the liquid and the gas colloidal phased. Strictly 
speaking, more sophisticated theories are needed in order to describe the interfacial properties near the triple point. 
DFT theories of the fluid-fluid interface in the two-component AO model suggest that for ideal polymers, the square 
gradient approximation underestimates the values of the interfacial tension 16 , and that it misses more subtle effects 
like oscillatory density profiles. Part of the difference come from using more accurate DFT's, but some also arises from 
the fact that the two-component DFT yields different phase-diagrams than the effective one-component description 1 ^. 
In particular, the two-component theories yield a larger distance between the critical point and the triple point 17,20 , 
which may have an important qualitative and quantitative effect on the interfacial behaviour. While it would clearly 
be desirable to have a two-component theory of similar accuracy for interacting polymers, this is not available at 
present. However, it is possible to use more sophisticated DFT approaches to study the one-component interacting 
polymer system^, a direction which we are now pursuing. In fact, we have already performed preliminary calculations 
using an accurate DFT for the HS part, and with the interaction of Eq. © treated as a perturbation. The results for 
surface tensions and interfacial widths are slightly higher than those from the present treatment, but the trends are 
very similar. This suggests that our use of the square gradient theory, coupled with our rater simple approximation 
for c(r), as used in Eq. JJjJ, is quite reliable. 

In the longer term, it would be interesting to develop some approximate density functional for interacting polymers 
using a two-component representation, perhaps along the lines of ref^£, in order to obtain accurate predictions of the 
interfacial properties of free fluid-fluid interfaces, adsorption and wetting behaviour at hard walls and even surface 
tensions and density profiles of fluid-solid interfaces. 

Although the actual values of the surface tensions found within the square gradient theory may only be accurate 
to about a factor of two, the differences between ideal and interacting polymers are fairly large, and follow from a 
simple physical explanation which seems robust. We therefore don't expect more sophisticated theories to reverse the 
trends discussed in this paper. On the other hand, whether or not the more subtle interfacial phenomena found for 
the two-component AO modeli 6 ^ will be even qualitatively similar for interacting polymers remains to be seen. For 
example, in the latter case the differences between the triple and critical points are less pronounced, which may lead 
to less well defined oscillations in the density profiles. Clearly more work remains to be done. 

The square gradient theory does not include the effects of capillary fluctuation*^, but these are not expected to 
be large, as shown in refA We can therefore make comparisons with experiments. Full two-component AO model 
calculations within DFT-t 6 " lead to values of the surface tension that are close to those of recent experiments 4 . However, 
our results suggest that including polymer-polymer interactions will lower the value of the surface tension, leading to 
less agreement. 



V. CONCLUSION 



We have used thermodynamic perturbation theory and the square-gradient approximation to calculate the prop- 
erties of the fluid-fluid interface for mixtures of colloids and interacting polymers within an effective one-component 
representation. We find significant differences compared to the case of ideal polymers. The main effect of polymer- 
polymer excluded volume interactions is to reduce the value of the interfacial tension 7 and the interface width W. 
This effect becomes more pronounced as the size-ratio q increases. It can be rationalised by the differences in depletion 
potentials: at the critical point, the range for the interacting polymer case is significantly less than for the ideal case. 
This has a pronounced effect on the square gradient prefactor k, and helps explain the differences between the two 
types of depletants. 
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FIG. 1: Phase diagrams of a colloid-polymer mixture for ideal (dashed lines) and interacting polymers (solid lines) compared 
for q = Rg/Rc = 0.67 (taken fromii). rj c is the colloidal packing fraction and rf v the reservoir polymer packing fraction. The 
binodal for interacting polymers is at a higher r\ v than that of ideal polymers, because the latter are stronger depletants than 
the former. 
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FIG. 2: Density profiles of the colloidal packing fraction for ideal and interacting polymers. In both cases, the size ratio is 
q — 1.05 and the relative polymer reservoir packing fraction from the critical point is (r/p — rf v ' cr%t } j rf v ' cr%t = 0.2. The polymer 
excluded volume interactions yield sharper interfaces. 
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FIG. 3: 10-90 width of the interface for ideal (white symbols) and interacting polymers (black symbols) versus the deviation 
of the polymer reservoir packing fraction from the critical point, for q — 0.34, 0.67 and 1.05. 
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FIG. 4: Dimensionless surface tensions for ideal (white symbols) and interacting polymers (black symbols) as a function of 
77c - r?f for q = 0.34, 0.67 and 1.05. 
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FIG. 5: Effective colloid-colloid depletion pair potential induced by interacting (solid line) and ideal polymers (dashed line) 
for q = 0.67 at their respective critical points, as a function of the distance between the particle surfaces x — r — 2R C . Inset: 
Ratio between the k parameter for interacting and ideal polymers at the critical point versus the size ratio q. 
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